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Abstract 



We prove that certain nonequilibrium expectation values in the boundary sine-Gordon 
model coincide with associated equilibrium-state expectation values in the systems which 
differ from the boundary sine- Gordon in that certain extra boundary degrees of freedom 
(^-oscillators) are added. Applications of this result to actual calculation of nonequilibrium 
characteristics of the boundary sine-Gordon model are also discussed. 
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1. Introduction 



In this paper we study tlie so called boundary sine-Gordon model with zero bulk mass 
(referred below as BSG). Its action is 

Absg = ^ dt dx{^l-^l) + - dtcos{^B + Vt) . (1.1) 

^'^9 J — oo J —oo 9 J — oo 

Here $ = $(x,t) is a scalar field defined on a half line — oo < x < 0, $b = $(0,t) is its 
boundary value, $t = 9t$, $x = and g, n and V are parameters. The way g enters 
( prrp allows one to interpret it as a quantum parameter, because it always appears in the 
combination gh; in what follows we set h = 1. In the quantum theory the parameter k 
carries the dimension S of [mass]^~^ . The boundary interaction in ( |1 . 1|) contains explicit 
time dependence through the term Vt. Although this time dependence can be eliminated 
by a simple change of the field variables — (which instead brings in the term 

linear in $t), the above form ( |1 . 1|) is more convenient for our analysis. Below we almost 
always assume that ^ > 0, more generally we will consider complex V with ^eV > 0. For 
future references let us write down the Hamiltonian corresponding to ( p. . 1|) , 



^^^^ 4ng 



1 /-o ^ 

— / dx{U^ + - -cos {^B+Vt) . (1.2) 
'^'^9 J —oo 9 

Here $(a:), H(x) are field operators obeying canonical commutation relations 

[U{x) , $(a;')] = -27rz£/ 5{x - x') (1.3) 

and again $b = $(0). 

At a nonzero V and a temperature T the system ( p..2|) develops a stationary nonequi- 
librium state which can be thought of as the result of an infinite time evolution of the 
equilibrium state of the corresponding "free" system, with the interaction term (the last 
term in ( |1.2[ )) adiabatically switched on. We will denote {A)bsg the expectation value 
of an observable A over this nonequilibrium stationary state. 

Besides being an interesting model of Quantum Field Theory on its own, the theory 
finds important applications in other branches of physics. As explained in , the model 
( pri| ) describes a quantum particle with the coordinate X = ^b + Vt and the potential 



^ We assume here that the normalization of the boundary field cos(^>_b + Vt) is fixed by the 
condition (cos($B(t) + Vt) cos(4>B(t') + Vt'))BSG 2"^ [i{t - t') + O)'^^ as < t - t' ^ 0. 
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energy —2nK cos{X), interacting with dissipative environment, the bulk part of the field 
$ playing the role of the latter (see also Iii this case V is interpreted as an external 

driving force. The model is also believed to describe an electric current through a 
point contact in the quantum Hall system ; in this case V is proportional to the voltage 
drop across the contact. In all cases the quantities of interest are the correlation functions 
of the boundary fields 

V+(t) = e^"^^(*) e'^' , V_(t) = e-'*«(*) e''^' . (1.4) 

Here $B(t) = S~^(t) $b S(t), and S{t) is the time evolution operator corresponding to the 
Hamiltonian (|1.2|). We will be particularly interested in the expectation values ( \± ) bsg- 
The main goal of this paper is to show that the expectation values ( \± ) bsg (and 
indeed some more general correlation functions of (|1 . 1| ) ) coincide with equilibrium expec- 



tation values of certain operators in a system which differs from (|1.2| ) in that it involves 
certain additional boundary degree of freedom. Namely, let us define a Hamiltonian 



H+ = 4^ c/x (H^ + $^ ) - y h - ^ (a_ e-- + a+ e"-- ) 



:i.5) 



where $(a;), Il{x) are again the Bose field operators obeying the same commutation rela- 
tions ( pr3|) as in ([L^), and the operators h, a_|_, a_ commute with $(a;), Il{x) and form 
among themselves the so called "g-oscillator algebra" , i.e. 

[h,a±]=±a±; q a+a._ - q'^ a_a+ = q - q'^ (1.6) 

with 

q = e^'^f . (1.7) 
Let be some representation of (|1.6D such that the spectrum of p+(h) is real and bounded 



from above. The Hamiltonian ( p..5| ) acts in the space 

n+ = j^0p+, (1.8) 

where JF is the space of states of the Bose field representing the commutation relations 
( p..3|) , and for V > this Hamiltonian is bounded from below. Then, for V > the system 
( |L5|) has a thermal equilibrium state described by the standard density matrix 

P+ = Z7i(K,y)e-™+ , (1.9) 
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where R is proportional to the inverse temperature, 



R = g/T 



and V) = Tr-^^ [e~'^^+ ] is the corresponding partition function. Here and below 

we treat g as a constant, and therefore we do not include it in the list of arguments of 
Let us denote (A)_(_ the expectation value of an observable A over this thermal 
equilibrium state. We will show that 



(V+)b5G = (W+)+ = (W_)+ 



where 



W+ 



w_ 



(1.10) 



(1.11) 



While the second equality in (|1.10| ) is a simple property of the equilibrium state ( pTOl ) 
(see Sect. 3), the relation between the nonequilibrium and equilibrium expectation values 
in (|1.10| ) looks rather unusual and suggestive. In fact, the Eq. (|l.lCl| ) is a particular case 
of a more general relation. Namely, the correlation functions in ( |1 . 2|) involving any num- 
ber of the Heisenberg operators V_|_(t) (but not V_(t)) coincide with the corresponding 



equilibrium-state correlation functions of the operators W_|_(t) in ( |1.5|) . 

Similarly, one can relate the correlation functions which involve the operators V_(t) 
(but not V_|_(t)) of ( pr2|) to certain equilibrium-state correlation functions. Let p_ be any 
representation of (|1.6| ) such that the spectrum of P-(h) is bounded from below. Consider 
the Hamiltonian 



H. 



a_ e 



) 



(1.12) 



acting in JF eg) p_ , and associated equilibrium state density matrix 



ZZ\k,V) e 



-i?H_ 



;i.i3) 



which is well defined for 3?e y > 0. Then 



(V_)bsg = (W_)_ = (W+). 



;i.i4) 



where ( A )_ stands for the expectation value of A over the equilibrium state (|1.13[ ), and 

a+e^^«; W_ = a_ 6"^^^ . (1.15) 
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Although the Hamiltonian (|1.12 ) is related to ( |1.5|) through a simple change V 



-V, we treat it as a distinct one because we always imply ^eV > 0. 

Besides being interesting by themselves, the relations ( p. .101 ) and ( |1.14| ) provide efficient 



tool for actual computation of the nonequilibrium expectation values {y^±) bsg- Using 
( p..lO|) , (|1.14| ) one can easily show that 



{V±)BSG=Td^logZ±{^,,V) , (1.16) 

where Z_|_ and Z_ are the partition functions in (|l]^) and ( |1.13| ) , respectively. On the other 
hand, the above partition functions are related in a simple way to the vacuum eigenvalues 
Q±(A,p) of the CFT analogs of Baxter's operators Q±(A) studied in (see Sect. 4 for 
some details), 

Z±{^,V)/Z±{0,V) = AT2-P//3'g±(A,p) . (1.17) 



where Z±(0,V) are the partition functions for the "free" systems defined by (|1.5| ) and 
( p..l2|) with K, = 0. The parameters A, p and (3 used in [^,0 are defined as 



A = m — '-^ (^] , p = -iV-^, (3=g. (1.18) 

Many exact results for these operators exist [^|^ . These include high and low temperature 
expansions, but most importantly, the vacuum eigenvalues Q±{X,p) are shown in to 
satisfy closed integral equations - the Destri-de Vega equations - which allow for their 
evaluation to an arbitrary degree of accuracy. In view of ( |1.10|) , ( [1.14| ) and ( |1.16| ) all these 



results directly apply to the non-equilibrium expectation values of the operators (|1.4|) . 
In applications, the quantity of a particular interest is 

J=V+{^,{0,t))BSG , (1-19) 

which is interpreted as the current through the point contact in the quantum Hall systemJl, 
or as the drift velocity in dissipative quantum mechanics. It was the main subject of interest 



^ The voltage V and current J ( |1.19| ) differs in normalization from the real voltage V^^^^'^^ 
and current J^^''^*) in the Hall system 

where e and h are the electron charge and Plank's constant. Also, g coincides with the fractional 
filling of the Luttinger state in a Hall bar and the temperature T is measured in energy units. 
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in many recent papers []5|,|8|,p|, p!0| , pl] ,|6| . In particular, for g = ^ in (|1.1|) , Bethe-Ansatz 
computation of this current was given in |p, under the assumption that the Boltzmann 
equation for the backscattering electrons in the Hall contact holds exactly. Further, it was 



conjectured in |11] that for all values of this parameter 



J = V + i7iTijdAog 



with 



/ g g 



g \27tT 



A-kT 



The "partition function" Z2p was defined in Ref. [|TT| as a power series 



;i.2o) 



(1.21) 



Z2p(^) = l + X]^'"/2n(p) 



;i.22) 



n=l 



with the coefficients l2n given by multiple infinite sums over all ordered sets m 
(mi, m2, . . . , rrin) of non-negative integers mi > m2 > ■ ■ ■ > m^ > 



1 



En 



r(mj + g{n — z + 1)) V{2p + m^ + g{n — z + 1)) 



^'"(^) m ^V[ + 9{n - ^) + 1) r(2p + m, + (/(n - + l) 



;i.23) 



Similarly looking conjecture 



J = V + i7rTKd^ log 



Z+{n, V) 
Z-{k, V) 



(1.24) 



where Z±{k, V) are just the partition functions in (|1.9| ), ( p..l3|) and ( |1.17| ), was proposed 
independently in 0. The (plausible) equivalence between ( p..20|) and (|1.24|) is not yet 
established. However, it is easy to see that ( |1.24D is a simple consequence of ( |1.16|) and 
the equation 

$,(0,t) = z7r«: (V+(t)-V_(t)) , (1.25) 
which is the boundary condition corresponding to (|1.1|), i.e. our result ( |1.1C1| ), ( |1.14| ) 



actually proves the conjecture (|1.24| ) 



2. Expectation values in boundary sine-Gordon model 

If y > the system ( p. . 1| ) evolves towards a stationary state which is characterized by 
nonzero expectation values ( $t ) bsg 7 ( ) bsg- This is not a thermodynamic equilibrium 
state and no simple explicit expression for its density matrix is known. In this section we 
describe the definition of this state in terms of real-time perturbation theory and discuss 
some properties of corresponding expectation values. 

Let us split the total Hamiltonian (|1.2|) into the free and interaction parts 



1 /-o ^ 

Ho=:^/ rfx(n2 + $2X cos{^B + Vt). (2.1) 

Using corresponding interaction representation one can write the density matrix P(t) as 

V{t) = e-^^o* S(t, -00) Po S(-oo,t) e^"°* , (2.2) 

where 

S(t,to) = Texpj dr Hf"*V)} = 

= 1+E(-^)' / ^^.({r})HS^'^*)(ri)HS^'^*)(r2)---Hf"*)(rO 

fc=l -^^o 

and 

HS*'^')(t) =e^"o* Hi e-'"o* . (2.4) 
In ( |2.3| ) and below the shorthand notation for the multiple ordered integrals 

T^kiir]) = I dn r dr2-- - r ' dTk (2.5) 

to -J to to J to 



2.3) 



is used. In writing (|2.2| ) we have assumed that the interaction has been adiabatically 
switched on in the infinite past (t = —00) when the system was in thermodynamic equi- 
librium at the temperature T, i.e. 

Po = Z^^ e-™o , R = g/T . (2.6) 

Then, the expectation value of an arbitrary operator A can be written as 

Tr^[P(t) A] = Tr^[Po S(-oo,t) A(*'^*)(t) S(t, -00) ] = ( S(-oo,t) A(*'^*)(t) S(t, -00) )o , 

(2.7) 



where (...)o denotes the expectation value over the equihbrium state (|2.6| ) of the free 
system, and the superscript "(znt)" means that this operator is taken in the interaction 
representation, i.e. A*^*"^*)(t) = g'^^o* A e~*^°*. Equivalently, one may write the above 
expectation value as 

(A)B5G = Tr^[PA(t)] , (2.8) 

where P stands for density matrix of the system at t = 0, i.e. P = P(0), and A(t) is the 
full Heisenberg operator 

A{t) = S(0, t) A^*'^*) (t) S{t, 0) . (2.9) 

All above formulae are very well known (see e.g. fl^); we included them here to fix 
the notations. In this paper we are interested in the expectation values ( V± ) bsg of the 
Heisenberg operators 

V±(t) = S(0, t) V5!"*)(t) S(t, 0) . (2.10) 

where 

V5!"*)(t) = exp{ ± z$g"'^(t) ± tVt} , (2.11) 

and $^"^*''(t) is the boundary field $5 in the interaction representation. It is convenient 
to introduce also auxiliary operators 

V±(t,to) = S(to,t)vf*^(t)S(t,to) , (2.12) 

where to is a parameter. For to = ( |2.12| ) coincide with the Heisenberg operators ( |2.10|) , 
and according to ( p.7| ) the expectation values of ( |2.10| ) can be expressed through ( |2.12| ) as 
follows 

(V,(t))BSG= lim (V,(t,to))o, (2.13) 

to — >— 00 

where a = ±1 . Our nearest goal is to prove the following useful representation for the 
operators ( ^.12| ) 

V.(t,to)=V(^"*)(t)x 

{l + E E C,(a|ai,...,a,) / D,({t}) V(^-*)(ti)---V(^f)(t,)}, ^ 

where the sum is taken over all arrangements of the "charges" ai, . . . ,ak = ±1, and the 
coefficients Cfc(a| ai, . . . , a/c) have the following explicit form 

k 

Cfc(cT|cTi,...,afc) = -j ai---akq ^3=1''^'^^ || sin(7r^?7j) , (2.15) 



i=i 



with 



i-i 

s=l 



(2.16) 



The fact that this operator can be written as the sum ( |2.14[ ) with some coefficients 
C'a;(c"| ci, ■ • • 7 cTfc) is obvious from its definition ( p.l2| ). Indeed, the series expansions ( |27^ ) 
for the evolution operators S in ( |2.12| ) allow one to represent ( |2.12| ) as a series of multiple 



integrals of certain products of . Then, using the commutation relations 

Vi'f\ti) Vi'^'Ht2) = V(^;*)(t2) ^^^T'\h) , ti > t2 , q = e^"" , (2.17) 

where (Ti,a2 — ±1, one can always rewrite each of these integrals as a combination of the 
time-ordered integrals as in (|T|). The easiest way to obtain the coefficients Ck in ( |2.14| ) 
is to notice that in view of ( p.3|) the operators ( p.l2|) satisfy the differential equation 

. d 



I — V^{t, to) = [ hS^"*) (to) , V<,(t, to) ] 

Oto 



with the initial condition 



(2.18) 



^2.19) 



It is easy to check that the expansion ( p.l4| ) satisfies ( p.l8| ) provided the coefficients Ck 
solve the recurrence relations 



Cfc(a| ai, . . . , (Tfc) = ^ (1 - Q-2-^^^) Ck-i{a\ ai, . . . , a^-i] 
^9 



(2.20) 



where the notation ( |2.16| ) is used. With the initial condition Co = 1, which follows from 
( p.l9| ), these relations lead to ( p.l5[ ). 

A simple consequence of ( 2.14 ) is the infinite series representation for the expectation 
value ( pT3D , 



2TiaT 



BSG 



K sin(7rgf) 



2n-l 



0"lv)f"2n-l J- 



J- A sm{TTg) 



(2.21) 

where the sum is taken over all arrangements of o"i, . . . , a2n-i = ±1 with zero total charge 
(J + Z]s=i ^ = 0, and 

Jn{(^o,cri, . . . , a2n-i\p) — 



n 

0<j<l<2n-l 



with To = 0. Here the parameters A and p are defined by ( |1.18| ). This representation (more 
precisely, the corresponding representation for the current ( |1.19D ) was previously obtained 
in Ref. by a combinatorial method. The above consideration provides an alternative 
derivation of this result. It also promotes the formula (|2.21|) to the level of the operator 
relation ( p.l4| ) which can be used in evaluating the multitime correlation functions. 

Let us emphasis here an important feature of the expression ( p.l4| ) and (|2.21|) . The 
sums over ai, a2, ■ ■ ■ <T2n-i there exclude configurations where any of the "cumulative 
charges" rjj, j = 1, . . . ,2n — l defined by ( |2.16| ) vanish, because one of the factors sin{7rgrjj) 
in (|2.21| ) (as well as in ( p.l5|) ) then turn to zero. As a result all integrals ( p.22|) appearing 
in ( |2.21| ) converge at the lower limit, and that is why the limit to ^ — oo in ( |2.13| ) poses 
no difficulty. And of course RHS of (|2.13|) , ( p.21|) does not actually depend on t as a 
consequence of the invariance of the expectation values upon an overall time shift. 



3. The system with g-oscillator 



Now we turn to the systems (T3) and ( |1.12| ) which involve the boundary g-oscillators, 
with the aim to derive our main relations ( p..lO|) and ( p.. 14 ). These two systems are quite 
similar and in what follows only ( |1.5|) is studied explicitly, and only (|1.10| ) is actually 
derived; the relation ( |1.14| ) can be obtained by obvious modifications of the arguments 
presented below. 

Like in (|T|), let us split (U) as H+ = Hq + Hi, 



Hn = Hr 



Vh. 



Hi 



K 



29 



(a_ e'"^^ + a+ e-'*^ ) 



(3.1) 



where Hq is given by ( |2.1| ) and the operators h, a± are defined in ( |1.6[ ). In the Matsubara 
representation the equilibrium density matrix (|1.9| ) reads 

(3.2) 



P+ = Z^\k, V) e-^"» S{-iR, 0) 

where the Matsubara operator S{—iR,0) is the "imaginary time" version of the time- 
evolution operator in the corresponding interaction representation, 



S(t2,ti) = Texp 



dt 



Hf"*)(t)} 



HS*"*^(t) = e^"°'Hie 



iUot 



-iUot 



Note that in view of the commutation relations 



H 



2/32 



(a_ V;!"*^(t)+a+VL'"*^(t)) 
9 



(3.3) 



(3.4) 



where V^"^*"* (t) are exactly the operators (|2.11| ) . 

We are interested in the expectation values of the operators (|1 . 1 1|) , which can be 
written as 

(W±)+ = Z+iTrH+[e-^"° S(-ii?,0) W±] . (3.5) 



The fact that ( W_|_)_|_ = (W_ )_|_ (as stated in ( 1.10|) ) is a particular manifestation of 



the detailed balance principle for the equilibrium system ( |1.5| ); on formal level it is easily 
established if one notices that the commutation relations are invariant with respect to the 
transformation a_|_ —>■ Aa_|_, a_ a_, where A is a constant. Therefore 

(W±)+ = Tc>, logZ+(/^,y) . (3.6) 

It is convenient for our purposes to introduce again an auxiliary time to and rewrite 
(P) as 

(W+)+ = Z+l'IVH+[e-^"°S(^o-^i^,^o)W+(0,^o)] , (3.7) 

where 

W+(t, to) = S(to, t) e^"o* W+ e-^«°* S(t, to) • (3.8) 



Our proof of (|1.10| ) will be based on remarkably simple representation for the operator 



W_|_(t, to) similar to ( p.l4| ). Using perturbative expansions for the operators S in ( |3.8|) with 



the explicit form ( |3.4| ) of the interaction Hamiltonian, and then applying the commutation 
relations ( |2.17| ) to achieve the full time ordering for the operators Vj!"^*\ one can bring 
( p78|) to the form analogous to ( p.l4|) , i.e. 

W+(t,to) = a_vf*^(t) + 

vf'^W E E C.(ai,...,a.) / P.({t})V(^;^*)(ti).--V(^f)(t.), 

k = l (Ti,...,crfc=±l -^^0 

where this time the coefficients Cfc(cri, . . . , a^) are not c-numbers, but some operators 
acting in In general, from this analysis one would expect these coefficients to be 
complicated polynomials in the operators a_|_ and a_. It turns out, however, that these 
coefficients contain only a_ but not a_|_, and moreover 

Ck{ai, . . . , cTfc) = (a_)^'=+i Ck{+\ ai, . . . , a^) , (3.10) 

where rjk is given by (|2.16|) with a = +1, and Cfc(a"| ai, . . . , a/j) are exactly the same 
numerical coefficients ( |2.15|) as in (|2.14| ). The rest of the notations in ( |3.9| ) is the same as 
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in ( 2.14 ). The derivation of ( |3.1C1| ) is very similar to our proof of ( 2.15 ) in the previous 
section. By its definition, the operator (|3.8|) must satisfy the differential equation 



i ^ W+ (t , to) = [ Hi'"*) (to) , W+ (t , to) ] 

Oto 



(3.11) 



with the initial condition 



W+(t,to)|,„=, = a_VV"'^^(t). (3.12) 
Substituting ( |3.9| ) with yet unknown operators Cfc, one obtains the recurrence relations 



I K 

2^ 



( Ca;_i(cti, . . . , crfc-i) a_ 



CA;-l(cri, . . . , CTfc-l) ) 

(3.13) 



with the initial condition Co = a_. It is straightforward to check that ( |3.1C1| ) satisfy 
( p.l3| ). Note that because of ( |3.10[ ) and (|2.15| ), the representation ( |3.9[ ) enjoys the same 
remarkable property as ( p.l4| ) , namely the series ( |3.9| ) contains only the terms where no of 
the "cumulative charges" rjj, j = 1, . . . , 2/c — 1 vanish. We want to stress also that this 
simple form of the coefficients Cfc is a very special property of the operator W_|_(t, to)- Had 
we taken, say, W_ instead of W+ in (P?^), no such simplification in ( |X^ ) would occur, 
i.e. the resulting coefficients Cfc would indeed come out to be complicated polynomials of 
both a_ and a+, and in particular the above "no zero cumulative charges" property would 
not hold. 



Our proof of (|1.10| ) is based on the representation (|3.9|) . Consider the expression 
for ( W+ ) and take the limit to — ^ — oo. Due to the above "no zero cumulative charges" 
property of (|3.9| ) the trace in (|3.7| ) in this limit factorizes as 



TrH+[e-™° S(to-zi?,to) W+(t,to)] 
TrH+[e-™o W+(t,to)] 



TrH+[e-™° S(to-ii?,to)] 



(3.14) 



and we obtain 



lim 

t(j — > — oo 



I^H+te-™" W+(0,to)] 



(3.15) 



The trace in the denominator is the partition function of ( |1.5| ) with k = 0. Now, using here 
the representation (|3.9| ) with (|3.10| ) for the operator W_|_(0,to), and taking advantage of 
the fact that 



e-™° (a_ 



(3.16) 
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one arrives exactly at the series in ( p.21|) for the expectation value (|3.15 ). This proves the 
relation (|1.10|) . Let us also note that using the same arguments based on the properties of 
( p.l4| ), (|3.9| ) and on (|3.16| ), one can establish more general relation between the multitime 
correlation functions of the systems ( p. .21) and (|1.5| ), namely 



;v+(ti) ■•• v+(t„))B5G = (w+(ti) -.■ w+(tO)- 



(3.17) 



Here W_(_ [t) = W_(_ (t, 0) are the Heisenberg operators associated with W 



+ ■ 



4. Discussion 

The relations ( |1.1C1| ) and ( |1.14| ) between nonequilibrium expectation values in ( |1 . Ij ) 
and equilibrium expectation values in ( |1.5| ) and ( |1.12| ) is the main result of this paper. 
We want to stress here that these relations are actually derived (in Sect. 2 and 3) from 
first principles, and in particular this result proves the validity of the expression ( |1.24| ) for 
the current ( p..l9|) conjectured earlier in p. General significance of this result depends on 
whether it can be extended to wider class of systems with nonequilibrium dynamics. At this 
time we feel reluctant to enter any speculations on this matter leaving that to when better 
understanding of the situation is achieved. Instead, in this section we discuss how these 
relations can be applied towards actual computation of the nonequilibrium expectations 
( V-i- )bsg and of the current (|1.19[) . 

According to ( |1.16|) and (|1.24|) , computation of these expectation values reduces to 



finding the partition functions Z±{k,V) of the equilibrium systems (|1.5|) and (|1.12[) . By 



itself this does not help much in solving the problem. Of course one can use Matsubara 
theory to write down a perturbative series in for Z±{k, V). However, n-th order coeffi- 
cients of this series are expressed in terms of 2n — 1 fold integrals (written down explicitly 
in [^) and it is not clear how to evaluate or simplify these integrals in general. In this 
respect this approach does not have any significant advantages over the series ( p.21|) ; in 
fact, ( |2.21|) looks more compact. The same remark seems to apply to the approach based 
on the representation ( |1.2C1|) conjectured in [11|. Even if ( |1.20 ) is proven, evaluation of all 



terms of the series ( |1.22| ) requires calculation of the sums ( |1.23|) of ever growing degree of 
complexity. Each of these representations — (|1.24| ) , (|2.21|) and (|1.20|) — provides more or 



less efficient way to calculate few first terms of the expansion of ( |1.19| ) but neither seem 
to give a full solution to the problem. In fact, if one is concerned with the actual evaluation 
of the coefficients of these expansions, even more efficient technique can be obtained 
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by combining certain analyticity properties and functional relations for Q operators (see 
below); we describe this technique in the Appendix. 

What makes the relations ( |1.1C1| ), ( p. .141 ) a useful computational tool is remarkable 
relation ( |1.17| ) between the partition functions Z±{k,V) and certain eigenvalues of so- 
called Q±-operators. These operators were introduced in ^ as CFT analogue of Baxter's 
Q-matrix |T^. The operators Q±(A) of ||] act in the Fock space of free Bose field with 



spatial coordinate compactified on a circle. The Fock vacua | p ) are parameterized by the 
value p of the zero-mode momentum. These vacua are eigenstates of the operators Q±(A), 
and we denote Q±{X,p) the corresponding eigenvalues, i.e. 

Q±(A,p) = (p|Q±(A)|p). (4.1) 

We will also use the notations 

A±iX,p) = X^''''^/^"Q±{X,p) . (4.2) 



The identity ([LTTD , i.e. 



Z±{k,V)/Z±{0,V)=A±{\p) , (4.3) 

(where the parameters k, V, g and A, (3"^ are related by ( |1.18| )) can be verified by direct 
comparison of the definition of Z±{k,V) in Sect.l and the definition of the Q-operators 
in [0. At the same time the above operators Q±(A) exhibit some very useful properties. 
First, the operators X'^'^'^'^^/f^ Q±(A), and hence the functions (^T|), are entire functions of 
A^, with known asymptotics at A^ — — oo along the real axis 0. Second, these operators 
obey certain functional equations (one is utilized in Appendix, see (A.5)), notably the 
famous Baxter's T-Q relation. These analytic characteristics and functional equations 
lead to closed nonlinear integral equation for the functions (|4.3|) (explicitly written down 
i), known as Destri-de Vega equation [1I^JT5[] . This equation has been used in to 



m 



analyze low-temperature expansions for the current (|1.19| ) and to confirm notable duality 
relation for this current first suggested in (some relevant results concerning this duality 
were obtained in |P,^|l6l). Detailed study of the solutions to this equation lies beyond the 
scope of this work. 



^ In the form written down in the Destri-de Vega equation applies only to the case of real 
2p > —g, i.e. pure imaginary V. For generic complex values of p the equation must be modified, 
namely the A integral over the real axis must be replaced by contour integral with the contour 
encircling all zeroes of A±{X,p) in the A plane, which become complex at complex p. 
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Appendix A. 

In this section we describe a technique of calculating the perturbative expansions in 
^2 (^"higji temperature expansions") of the expectation values ( p..lO| ), ( |1.14| ) and ( |1.19| ). 
The technique is based on the functional relations for the Q-operators derived in 1^,0. 
Unlike the main text of the paper which was meant to be more or less self-contained, this 
Appendix relies heavily on the results of Refs. , so some familiarity with these works 
will be useful in reading it. 

Let us remind some properties of the functions ( [4.2|) which follow from the analysis 



(i) The functions A±{X,p) are entire functions of A^, and A±{0,p) = 1. The power 
series @ , 



n=l 

has finite radius of convergence which is determined by position of closest to origin zero 
of A_|_(A, p). For real 2p > —g this zero is a real and positive. Let us mention here that 
according to (|1.16| ), ( |2.21| ) and (^73|), the coefficient an{p) can be expressed in terms of the 



in 



oo 




(A.l) 



integrals (|]2|), 




2n-l 



The coefficients an{p) here differ in normahzation from (p) used in 



[Q, namely 



(p)=5-'" T'-il-g) hL-''^\p) . 
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with CTg = ±1, ry-, = 1 + X]^=i "^s- Note that this relation holds as is only when < 
g < ^ and 3fJe (2p) > —g, since otherwise the integrals ( p.22| ) generally diverge; analytic 
continuation in these parameters is required outside this domain. 

(ii) The coefficients an{p) are meromorphic functions of p. Moreover, they are analytic 
in the half plane 3f?e (2p) > —g. The last statement follows directly from Eqs. (|2.22|) and 

(iii) The following asymptotics, 

aAp) ^ f ' "'^ r^"(l - 9) ■ a. p^oo, (A.3) 

hold in the half-plane 3?e (2p) > —g. 

(iv) A_(A,p) is related to A+(A,p) as 

oo 

A_(A,p) - A+(A, -p) = expj - <^n{-p) A^" } . (A.4) 

n=l 

(v) The functions A±{X,p) obey so called "quantum Wronskian" condition, 

^2nip A+{qH,p) A_{q-^X,p) - e-2^*P A+{q-^X,p) A_{qH,p) = 2i sm{2np) , (A.5) 
where q = e^'^^ . 

The Eqs.( |A.4D and ( [A.5| ), supplemented with the analyticity (ii) and asymptotic con- 
ditions (|A.3| ) constitute a Riemann-Hilbert problem which defines the functions A±{\,p) 
completely. Indeed, substituting the expansions ( |A.1|) , (|A.4|) into the ( |A.5|) , one obtains 
relations of the form 

sin{Txng + 2'np) a^ip) - sm.{Tmg - 2Txp) an{~p) = Rn{p) , n = 1, 2, . . . , (A. 6) 

where the functions Rn{p) are expressed through ak{p) with k = 1, . . .n — 1 only. For 
example, 

i?i(p) = 0, 

R2{p) = -{qai{p) + q-^ai{~p)Y e^'^^ sin(27rp)/2 . 
Since an{p) are analytic in the half plane 3?e {2p) > —g, and an{p) const pi-2n+2ng 
p — s> oo there, one can solve the Eq.( |A.6|) with respect to an{p), 
nT{l - 2g) V{g + 2p) 



ai{p) 



sin{Tig) r(l - 5f + 2p) 



/ ^ / -,\n ^ T{2-n + ng ■^2p) 

aM = -ir - TT. TTTX A.^ 

n L [n — 1 — ng + 2p) 



+ 00 



dx Rniix) V{n — 1 — ng + 2ix) r{n — 1 — ng — 2ix) 
2n X + ip 
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where n = 2, 3, . . . .In writing ( |A.8|) , we assume that < g < ^ and 3?ep > 0. As the 
functions Rn{p) are uniquely expressed (by the use of ( [A.7| )) through "lower" ak{p) with 
/c = 1, 2, . . . , n — 1, the Eq.( |A.8|) provides a recursion for evaluation of an{p)- It allows one 
to represent an{p) in terms of {n — l)-fold integral. Let us present explicit formulae for 02 

and as {0 < g < ^ and 3f?ep > 0), 



dx Si(x) 



- 2 

^3 



2% X + ip 



^ ^ r3(i+^7) r(2-3(7 + 2p) ^ ^ 

sin(47r^) /•+°° rfy /'+°° rfx S2ix,y) ^1 [^^ dx Ss{x) 



7r2 271- 27r (y + - y - zO) 3 i-oo 27r x + 

where the functions S are 

^i(x) = sinh(27rx) r(l -2g + 2ix) r(l - 2^? - 2ix) {r{g + 2ix) r{g - 2ix)f , 

^2(0;, y) = sinh(27r|/) sinh(27rx) r(£? + 2iy) r{g - 2iy) r{2g + 2iy) r{2g - 2iy) x 

r(2 - 3^ + 2iy) r(2 - 3^ - 2iy) V{1 ~ 2g + 2ix) V{1 - 2g - 2ix) x 

{T{g + 2ix)T{g-2ix)f , ^^'^^^ 

Ss{x) = sinh(27rx) r(2 - 3g + 2ix) r(2 - 3g - 2ix) {r{g + 2ix) T{g - 2ix)f x 
sin(4i7rx + 2'n:g) — 2 sin(27r(7) 
sin(2z7rx + 2'n:g) 

Eqs. (|A.(j| ), ( |A.y| ) have obvious computational advantages over ( [A.2| ). First, the number 
of integrations is greatly reduced - while (|A.2| ) requires evaluation of 2n — 1-fold integral 
( p.22|) to compute an{p), the recursion relation ( |A.8|) leads to n — 1 fold integral for this 
quantity. Second, it is very convenient in numerical calculations because the integral in 
(|A.8|) converges very fast at infinity. It is important to note also that analytic continuation 



of the expressions (|A.8| ), (|A.9D outside the domain < g < ^ and 3fJep > is rather 
straightforward — it is done by appropriate deformation of integration contours. For 
example, the analytical continuation of 02 (p) into the domain ^ < g < 1 and Sftep > 
reads, 

a M ^2^-^9 r(2^ + 2p) r /•+°° dx 

' r2(i + g) r(l -2g + 2p) 1 27, x + ip 



2 ' -V- i w-00 '271 x + ip 

sm{27ig)T{3-4g)T^{l-g) T^{3g 



{2p + l-2g) {2p-l + 2g) 
where the function Si{x) is the same as in ( |A.10| ). 
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